Ultracold Li + Li2 collisions: bosonic and fermionic cases 
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We have carried out quantum dynamical calculations of vibrational quenching in Li + Li2 collisions 
for both bosonic and fermionic ''Li. These are the first ever such calculations involving fermionic 
atoms. We find that for the low initial vibrational states considered here [v < 3), the quenching 
rates are noi suppressed for fermionic atoms. This contrasts with the situation found experimentally 
for molecules formed via Feshbach resonances in very high vibrational states. 
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There is at present great interest in the properties 
of cold and ultracold molecules. At the end of 2003, 
three groups Q, |^ |^ created long-lived molecular Bose- 
Einstein condensates (BECs) from ultracold atomic gases 
by making use of magnetically tunable Feshbach reso- 
nances. The molecules are formed in highly excited vibra- 
tional states, and any vibrational relaxation will result in 
trap loss. A crucial breakthrough was the use of fermionic 
isotopes of the atomic gases (^Li, ""^K) tuned to large pos- 
itive scattering lengths; it was found that in this case the 
atom-molecule and molecule-molecule inelastic collisions 
are strongly suppressed 0, H, S 0i and the molecular 
cloud consisting of weakly bound dimers of fermions ex- 
hibits a remarkable stability against coUisional decay. By 
contrast, for weakly bound dimers of bosonic atoms the 
vibrational relaxation is fast, so that the corresponding 
molecular cloud decays quickly and quantum deg eneracy 
can be achieved only transiently |^ l^, 0, 0, Petrov 
et al. |l3l| have explained the important difference be- 
tween molecular gases of weakly bound dimers formed 
from bosonic and fermionic atoms in terms of different 
quantum statistics for the atoms. 

In this Letter, we carry out the first full quantum 
dynamics calculations for spin-polarized collisions of ul- 
tracold homonuclear molecules for Li + Li2 in both the 
bosonic (^Li) and fermionic (^Li) cases. We find that for 
low-lying vibrational states there is no systematic sup- 
pression of the vibrational relaxation (quenching) rates 
for fermions, even when the scattering length is large and 
positive. The analysis of Petrov et al. 13] applies only 
when the scattering length is large and positive and the 
molecules are in their highest-lying vibrational state. 

Quantum scattering calculations require accurate po- 
tential energy surfaces. We have previously shown that 
the interaction potential for spin-polarized (quartet) Li + 
Li2 is highly non-additive [15| , with a well depth 4 times 
that of the sum of Li-Li pair potentials. In the present 
work, we use a global Lis potential obtained from all- 
electron coupled-cluster electronic structure calculations, 




FIG. 1: Cut through Lis quartet surface at coUinear geome- 
tries, showing seam between and ''ll states. Contours are 
labelled in cm~^. The two minima are at depths of 760 and 
950 cm-^ 



which will be described in detail elsewhere [l^ . In par- 
allel work, Colavecchia et al. ^] have obtained another 
potential surface based on pseudopotential calculations. 
However, we chose to use our own surface in the present 
work because it takes better account of some features 
of the potential. In particular, as pointed out in Ref.Hsl 
there is a seam of conical intersections at collinear geome- 
tries, between the surface that is the ground state at 
long range and a *n surface that becomes the ground 
state inside ri = r2 = 3.1 A. The resulting cusp drops to 
energies just below the atom-diatom threshold, as shown 
in Fig. n so may play a significant role in the dynamics. 
Our potential takes better account of this seam than that 
of ref . 

The methods we use to carry out quantum dynamics 
calculations on systems of this type have been described 
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in our previous work on Na + Na2 0, and elsewhere 
in the context of thermal reactive scattering so a 
brief summary will suffice here. The potential energy sur- 
face is barrierless, so that it is essential to take reactive 
(atom exchange) collisions into account. The positions 
of the nuclei are described in hyperspherical democratic 
coordinates. The configuration space is divided into in- 
ner and outer regions, and the boundary between them is 
placed at a distance such that couplings due to the resid- 
ual atom-diatom interaction can be neglected outside the 
boundary. In the inner region (hyperradius p < 45 ao in 
the present case) , we obtain the wavefunction for nuclear 
motion by propagating a set of coupled equations in a 
diabatic-by-sector basis that is obtained by diagonalizing 
a fixed-p reference Hamiltonian in a basis set of pseudo- 
hyperspherical harmonics. In the outer region, we use 
the Arthurs-Dalgarno formalism |21(|, which is based on 
Jacobi coordinates, and compute by inwards integration 
regular and irregular solutions of a radial Schrodinger 
equation which includes the isotropic {R^^) part of the 
interaction. Matching between wavefunctions in the in- 
ner and outer regions yields the scattering S-matrix. 

The difference between the ^Li and ^Li cases is imple- 
mented by selecting the pseudohyperspherical harmonic 
basis functions included in the basis set. For both bosons 
and fermions, we consider spin-stretched states of all 
three atoms involved (states with \Mf\ = F = i^max)- 
For such states the nuclear spin wavefunction is sym- 
metric with respect to any exchange of nuclei. However, 
the electronic wavefunction for the quartet state is anti- 
symmetric with respect to exchange of nuclei. Thus, to 
satisfy the Pauli Principle, the wavefunction for nuclear 
motion must be antisymmetric with respect to exchange 
of bosonic nuclei (for fermionic alkali metal atoms) and 
symmetric with respect to exchange of fermionic nuclei 
(for bosonic alkali metal atoms). 

Asymptotically, the hyperspherical functions correlate 
with atom-diatom functions; the diatom functions for 
Hund's case (b), which is appropriate for L\2 (^Sj), are 
labelled with a vibrational quantum number v and a me- 
chanical rotational quantum number n; n couples with 
the diatomic electron spin s = 1 to give a resultant j, but 
for Li2 the splittings between states of the same n but 
different j are very small and are neglected here. In the 
state, only even n is allowed for ^Li2 and only odd 
n for ^Li2. In the outer region, we include rovibrational 
states with v = 0, 1, ...,7 with all even rotational levels 
up to 7imax = 32, 30, 28, 24, 22, 18, 14, 10 for bosonic ^Li2 
and rovibrational states with v = 0, 1, 7 with all odd 
rotational levels up to rimax — 31, 27, 25, 23, 19, 17, 13, 7 
for fermionic ®Li2. In the inner region, the number of 
coupled equations for bosonic atoms varies from 97 for 
total angular momentum J = to 827 for J = 10 and 
for fermionic atoms it varies from 85 for J = 0+ to 782 
for J = 11". 

We have calculated the low-energy elastic and total 



quenching cross sections a-{E) for ''Li2 colliding with ^Li 
and for ^Li2 colliding with ^Li, for bosons and fermions 
in all allowed initial states (u, n) with v < 3 and n < 11. 
The cross sections for v = 1 (with n = for ^Li2 and 
rt = 1 for ^Li2) are shown as a function of collision energy 
E in Figs. 121 and 131 together with the contributions from 
individual partial waves J. The partial wave sums are 
converged up to 500 mK for Jmax = 10. For v = 1, both 
the elastic and inelastic cross sections are slightly smaller 
for ^Li than for ^Li at energies below 1 mK. However, the 
difference is only about a factor of 2, and it is not due 
to the difference between boson and fermion symmetries: 
for other values of v and n the ratio is often reversed. 

At very low energies (up to about 1 mK), the cross 
sections are dominated by collisions with orbital angu- 
lar momentum I — 0, for which there are no centrifugal 
barriers. A consequence of the symmetry constraints is 
that total angular momentum J = has an ^ = chan- 
nel for ^Li -I- ^Li2 but not for ^Li -I- ^Li2; in the latter 
case the lowest partial wave containing ? = is J = 1~, 
where the superscript indicates odd parity. For this rea- 
son, low-energy calculations on ^Li are significantly more 
expensive than those on ^Li. 

For nearly all initial v and n, and for both bosons and 
fermions, the ratio of inelastic to elastic cross sections is 
greater than 30 at i? = 1 pK and remains greater than 1 
throughout the Wigner regime (i.e., up to about 1 mK). 
The only exceptions to this are for rotational quenching 
of the {v,n) = (0, 2) state for bosons (ratio 4.3 at 1 /iK) 
and the (0,3) state for fermions (ratio 7.5), for which 
there is only one inelastic channel with ^ = 0. 

Scattering calculations for Li -I- Li2 are converged up 
to about lO"'^ K using just a single partial wave ( J = 
for ''Li and J = 1~ for ^Li). At higher energies, other 
partial waves start to contribute significantly. For ^Li 
the contributions decrease monotonically with J at low 
energies because of centrifugal barriers corresponding to 
Z = J torn = 0. For ^Li, the situation is slightly different: 
each partial wave is made up of a "parity favored" block 
(with parity (—1)-') and a "parity unfavored" block; the 
former includes an n = 1 channel with I — J — 1, while 
for the latter the lowest I is I — J. 

At collision energies below the centrifugal barrier, the 
partial cross sections follow Wigner laws, crl^i^i ^ E^~^^'^. 
Since hnci{E) = i;coiicri„ei(-B), where Vcoii = {2E/^Y/'^ is 
the collision velocity and fx is the atom-diatom reduced 
mass, the inelastic rate coefficient is independent of E 
at limitingly low energies. For initial v = \, the results 
in Figs. 121 and give calculated limiting values fcind = 
5.6 X 10^^° cm-^ s~^ for bosons and 2.8 x 10~^" cm'^ s~^ 
for fermions. 

At energies above the barrier height, the inelastic prob- 
abilities saturate at values close to 1 and the energy de- 
pendence in each partial wave is then governed by the 
prefactor E~^ that multiplies the S-matrix term in the 
expression for the cross section. At energies high enough 
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FIG. 2: Elastic cross sections (upper panel) and inelastic cross 
sections (lower panel) for ^Li + ^Li2(u,; = 1, rii = 0), with 
contributions from individual partial waves. The vertical lines 
indicate centrifugal barrier heights for Z > 1. 



FIG. 3: Elastic cross sections (upper panel) and inelastic cross 
sections (lower panel) for ®Li + ®Li2(wi = 1, rii = 1), with 
contributions from individual partial waves. The vertical lines 
indicate centrifugal barrier heights for / > 1. 



that several partial waves (typically > 3) contribute sig- 
nificantly to the overall cross sections, the energy depen- 
dence is described approximately by classical Langevin 
capture theory j2^ . 

hUE)^M'^E/^,Y'^(^^'^'\ (1) 

where Cg = 3086 -E/iOq is the isotropic atom-diatom dis- 
persion coefficient. The Langevin result is compared with 
the full quantum results for the quenching rates for v = 1 
and 2 in Fig. 0] it may be seen that the Langevin and 
full quantum results agree semi-quantitatively at colli- 
sion energies above about 10 mK for both bosons and 
fermions. 

The atom-atom part of the potential used in our cal- 



culations gives triplet scattering lengths that are approx- 
imately correct for both ''Li (—26.35 ag) and ^Li (—1531 
oq). However, at the magnetic field used in the molecule 
production experiments 0], the spin states involved had 
a w -f3500 ao, and it is important to establish whether 
fermionic suppression of quenching rates occurs for such 
large positive atom-atom scattering lengths. We have 
therefore repeated the calculations for energies E = 1 /iK 
and 1 nK with the short-range part of the 2-body poten- 
tial adjusted to reproduce a ~ 3500 oq. (Slightly different 
adjustments were needed for bosons and fermions). Un- 
der these circumstances the calculated quenching rates 
for initial w = 1 at limitingly low energy are 7.7 x 10^^^ 
cm^ s-i for ^Li and 2.9 x 10"!° cm^ s'^ for ^Li. This 
confirms that there is no fermionic suppression for the 
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FIG. 4: Total inelastic rate coefficients for collisions of Li 
with Li2 {v = 1 and 2, with n — for bosons and n — 1 for 
fermions) . 



low-lying states that we consider. 

The Langevin result, Eq. ^ may be used to predict 
inelastic rates for other atom-diatom alkali systems in 
spin-polarized states. Outside the Wigner regime, the 
inelastic rate coefficients for different systems will be pro- 
portional to I [1^1''^ . The resulting inelastic rate co- 
efficients for ^^Na, '"'K, ^'''Rb and ^^^Cs are lower than 
that for lithium by factors of 1.81, 1.75, 2.43, and 2.65, re- 
spectively, using Cg coefficients equal to twice the atomic 
dispersion coefficients from Ref. |23. The lower bound of 
applicability of the model may be estimated from the 
centrifugal barrier heights for Z = 3, which are 6.86 niK, 
1.89 mK, 0.537 mK, and 0.235 mK, for Na, K, Rb, Cs, 
respectively. 

In conclusion, we have carried out full quantum dy- 
namics calculations of vibrational relaxation (quenching) 
of ''Li2 by ^Li and of ^Li2 by ^Li, taking full account of 
the boson and fermion symmetries. We find that for low 
initial vibrational states of the molecules [y < 3), there 
is no systematic suppression of the quenching rates for 
molecules formed from fermionic atoms, even when the 
atom-atom scattering length is large and positive. This 
contrasts with the situation for molecules formed in the 
highest vibrationally excited state Q, |^ |^ . 

Collisions of molecules involving mixed isotopes of- 



fer fascinating new possibilities and in some cases allow 
identifiable chemical reactions even between ground-state 
molecules. For example, the reaction ^Li -I- ®Li2(w = 0) 
^Li -I- ^Li^Li(i; = 0) is exothermic by 2.643 K because 
of the difference in zero-point energy of the two diatomic 
molecules. Calculations on collisions such as this will be 
described in future work. 
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